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The NLO generalization of the dipole factorization formula for the structure functions F2 and 
Fl at low X is calculated using light front perturbation theory. That result gives some interesting 
insight into the kinematics of initial state parton showers in mixed space. 



I. INTRODUCTION 



QCD at low Bjorken x has become a very active topic of research in particular thanks to deep inelastic scattering 
(DIS) experiments at HERA. A large part of the DIS phenomenology at low x since the start of HERA has been 
based on the dipole factorization derived by Nikolaev and Zakharov in Ref.[l|. Contrary to earlier works in low x 
O I QCD, which have been performed in momentum space like the derivation of the BFKL evolution 043 resummming 
I . the high-energy leading logs (LL), the dipole factorization is formulated in mixed space, specifying the transverse 
position and the light-cone momentum of the involved particles. The dipole factorization has been proposed not onl 
(— I I for DIS structure functions but also for other DIS observables at low x, like diffractive structure functions \m 
I— I. deeply virtual Compton scattering and exclusive vector meson production Q. The dipole factorization provides a 
very intuitive picture of DIS observables: the virtual photon radiated from the lepton first fluctuates into a quark- 
, antiquark dipole, which then interacts with the target via gluon exchange(s). In the leading order (LO) version of 
^ ^ • the dipole factorization, the splitting probability for the photon into a dipole is known from perturbative QED. By 
[ contrast, the other factor, which is the dipole-target elastic scattering amplitude, contains all the QCD dynamics, 

■ both perturbative and non-perturbative. 

■ Numerous phenomenological studies based on the dipole factorization have been performed. In most of them, 
. phenomenological models for the dipole-target amplitude encoding various effects have been used. However, it has 

been soon realized that the BFKL evolution can be rederived as the low x evolution of a dipole cascade in mixed space 
0, S]. Hence it is very natural to combine the two results, and pick the dipole-target amplitude among the solutions 
of the BFKL equation in mixed space. By adding such constraint from perturbative QCD, one reduces a priori the 
needed amount of modeling, down to the choice of the initial condition for the BFKL evolution. 

However, the BFKL evolution has some severe shortcomings like the violation of unitarity at high energy and the 

■ sensitivity to the non-perturbative infrared physics, especially if the coupling is running. The phenomenon of gluon 
saturation at high energy [1, is both a consequence of those issues and a quasi-perfect solution^ to them. When 
taking gluon saturation into account, the BFKL equation is generalized into the B-JIMWLK equations, derived both 
from the high-energy operat or p roduct expansion of Wilson line operators and from the Color Glass Condensate 
(CGC) effective theory (l5l - l22j based on earlier works [23l - [25| . In a mean-field approximation, the B-JIMWLK 
equations reduce to the Balitsky-Kovchegov (BK) equation, also derived [H, [131 independently in the framework of 
Mueller's dipole cascade It has been soon realized that when using the LO dipole factorization together with the 
LL BK or B-JIMWLK equations, one describes the DIS data only qualitatively, because the obtained low x evolution 
is faster than in the data. One is thus led to consider higher order corrections. 

As a first step towards the NLO/NLL accuracy, the first contributions to the running of the coupling as in the 
BK or B-JIMWLK equations have been calculated [28'-'3^, leading to appropriate prescriptions to set the scale of the 
running coupling in the BK and B-JIMWLK equations. By simply promoting the coupling the BK and B-JIMWLK 
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^ Gluon saturation solves the problem of infrared sensitivity at large enough rapidity and restores the unitarity of the dipole-target 
amplitude at fixed impact parameter. However, it is not enough jlll. Il2ll to prevent the violation of the Froissart bound [T3II . after 
integration over the impact parameter. 
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equations to a running coupling following those prescriptions, an very good description of the low x DIS data can 
already be achieved j3ll - [33| . 

Then, the full NLL BK equation has been calculated [s^, UHl- Finally, the NLO photon impact factor has been 
obtained [SG^ in full coordinate space in the high-energy operator expansion of the product of two electromagnetic 
currents. That coordinate space representation is especially convenient to study conformal symmetry breaking or 
restoration depending on the choice of factorization scheme for the LL resummation. That NLO calculation allows to 
check explicitly that the B-JIMWLK equations are properly resumming the LL. Notice that the ongoing effort towards 
NLO accuracy with gluon saturation effects is not limited to DIS observables. As an example, NLO corrections to 
inclusive forward hadron production in hadronic collisions has been calculated most recently [37| . That observable is 
probably the simplest one for hadronic collisions to be sensitive to gluon saturation. It also involves an effective gluon 
distribution directly related to the same dipole target amplitude as used for DIS observables. 

Unfortunately, however, the results [s^ for the NLO photon impact factor are not available in a form suitable 
for phenomenological studies. One needs to perform non-trivial Fourier transforms from full coordinate space to full 
momentum space for the incoming and outgoing photons in order to obtain the NLO generalization of the dipole 
factorization formula I*]. The purpose of the present paper is to fill that gap. More exactly, the explicit Fourier 
transform of the previous results [3g is left for further studies, and instead the NLO generalization (|58p of the dipole 
factorization formula for the virtual photon cross sections a^{x,Q'^) and aJ^{x^Q^) is directly calculated. Let us 
remind that those cross sections are related to the DIS structure functions Fl and F2 = Ft + Fl as 

The calculation presented here is based on the same general ideas as in ref. [SG^ , however the method differs significantly. 
Indeed, light-front perturbation theory 38, 39] is used here whereas the more standard covariant perturbation theory 
is used in ref. [36j . 

The results of both ref.fs^ and the present paper should give, in the two-gluon exchange approximation, the 
appropriate NLO photon impact factor for the BFKL equation in mixed space^ The NLO photon impact factor for 
the BFKL equation in momentum space has also been calculated previously (4Cll442l | . but is not available in a closed 
analytical form. The comparison of the three approaches is left for further studies. 

In the section |lll the calculation of the quark-antiquark-gluon components of the light-front wave functions of 
transverse or longitudinal virtual photons is performed, leading to the results ([T2|) . ((T3)) and ([T4)). The associated 
virtual corrections are inferred from probability conservation. Those results are buildings blocks not only for the NLO 
corrections to the structure functions, but also to other DIS observables admitting a dipole factorization. 

Those intermediate results also allows to shed some light on the kinematics of dipole cascades in mixed space. A 
simple prescription (j24p to take recoil effects into account directly in mixed space is proposed. And a very simple 
expression in mixed space (j29p for the formation time of arbitrary multi-parton Fock states through initial state 
radiation is conjectured. That expression is purely kinematical and depends on the result of the parton cascade but 
not on the diagram followed. That expression should nevertheless contain full recoil effects diagram by diagram. 

The results of the section |lT] are used in the section IIIII to derive the NLO generalization ([55)) of the dipole fac- 
torization formula for the virtual photon cross sections a]^{x,Q^) and a^{x,Q'^). That expression involves the NLO 
impact factors ([501) (l5Tt . Conclusions are then given in the section ITVl 

For completeness, a detailed description of the formalism and notations used for the calculation is provided in the 
appendix [Xj 

II. VIRTUAL PHOTON WAVE-FUNCTIONS 

As explained in appendix the first step in the calculation of the DIS cross section is to obtain the light-front 
wave-functions of longitudinal and transverse virtual photons at the appropriate order. Notice that if the photon 
stays a photon or splits into a lepton pair, it cannot interact with gluons of the target and thus does not give a 
large contribution to the DIS cross section at low x. So, we will always consider only the components of the photon 
wave-functions containing colored partons, and drop the other ones. 

A. Quark-antiquark components at LO 

As a warm-up, let us outline how to obtain the standard LO results for the virtual photon wave-functions in the 
formalism based on light- front perturbation theory [H, [s^ described in appendix [Xj At this order, the only diagram 
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FIG. 1: LO diagram for the transverse photon wave function. 



contributing to the colored sector of the transverse photon wave-function is shown on FiglTJ For a photon of hehcity 
A = ±1, virtuahty Q, momentum large enough and q = one gets from the general formula (jAlOp the incoming 
state in the Heisenberg picture 

|7t('Z+,0',A)^)^o = 5]fe^(ko,fco+,/io,Ao,/)dt(ki,fc+,/ii,Ai,/)|0)-^^^ \ — 

states { ~ 2q+ ^0 ^1 1 

X (0|d(ki, /ii, Ai, /) &(ko, /lo, ^0, /) Um ^7^(0, A)|0) 

'■^ dzo dzi \ r fnu f • ko 



Jo V^O Jo v^i . J 



2(2^)70 y^Jo V^i ' ' 'J (2^)2(zoziQ2+ko') 

X E E ^^e^O' ^0+' Ao, /)dt(-ko, k+, -ho, Ao, /)|0) , (2) 

/ Ao 

where Zi = /q'^ is the longitudinal momentum fraction with respect to the momentum of the virtual photon. 

For the case of longitudinal photon, the LO diagram for the wave-function is essentially the same as on Figdl but 
with the effective vertex (jA2ip . so that 

|72(9+, = E ^^(^0, ■ ■ ■ )d^(ki, . . . )|0) ^ V^j(,+ ,Q.)^,(k„....) ,(k„...) 



qq states ^ 2q+ ^0 ^1 

- e dzp dzi f (fko zpZiQ 

{27r)i ^Jo V^i ' 'J {2nrizoz,Q^ + W) 

^E^/EE b^(^o,k+, ho, Ao, f)d\~ko,k^ ,—ho, Ao, f)\0) . (3) 

/ Ao ho 

Performing Fourier transforms (|Alip thanks to the relations (|Bip and (|B2I) , one obtains the mixed space represen- 
tation of these wave-functions, which write 

ho 

^E^/E bH^o,k+,ho,Ao,f)dH^i,k+,-ho,Ao,f)\0), (4) 

/ Ao 

with 



$r°(xo,xi,zo,zi,/io, = i[zi-zo- {2ho)X\ Q\l zozixl^ Ki ( zozixl^ ) (5) 



i) = -2zo ^^1 Q Ko (^\J zozixl^ , 



^i*^ xo,xi,zo,zi = -2zoziQKo\Q\J zozixl^] , (6) 



using the notations Xpg = Xp — Xq and = |xpq|. And Ko(a;) and Ki(a;) are modified Bessel functions of the second 
kind. Those LO results are of course not new, and are consistent with previous calculations in the literature like in 
ref.El. 
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-< ki or xi, hi, Ai, f 




q+.Q.X 



^ k() or X|-,, A:(J", h^., A,,, / 



or X2, . A2, a 



ki or xi, k^, hi, Ai, / 



^AAAAAAAAA> 



k|, or xq, k^., hi,.. Aq, / 

ko or X2. A;^, A2, a 




ki or xi, A;^, /ii, Ai, / 



^AAAAAAAAA4 




k|, or X||, fc^, /lo, A|), / 



k^ or X2, k^ , A2, a 



ki or xi, fc^, hi, Ai, / 



FIG. 2; Diagrams for the quark-antiquark-gluon component the transverse photon wave function at NLO. Each diagram is 
ordered along x'^ , from — s- — oo on the left to = on the right. 



B. Quark-antiquark-gluon components in the NLO wave-functions 



At NLO, the colored sector of the virtual photon wave- functions contains both quark-antiquark and quark-antiquark- 
gluon components, so that on can write 



/ NLO ' QQ ' ' 



qqg 



(7) 



The four light- front diagrams on Figl^^a), ^b), [5Jc) andj^d) contribute to the quark-antiquark-gluon component 
of the transverse photon wave- function at NLO. In the diagrams [D^ a) and^b), the gluon is emitted from the quark 
or from the anti-quark at a later than the splitting of the photon into the quark-antiquark dipole. By contrast, 
in the diagrams Wis) and EJd) , the splitting of the photon into a quark-antiquark-gluon system is instantaneous in 
a;"*". Those two different instantaneous interactions correspond to the two terms in the instantaneous vertex (lAlSp . 
Hence, one has 



|7;(9+,Q',A)^)^ 



E 



6t(ko,...)dt(ki,...)at(k2,...)|0) 



qqg states 



E 



2q+ '^0 '^1 '^2 + 

(0|a(k2, . . . )b{ko, . . . )Ui{0)b^{k,, . . . )|0) (0|d(ki, . . . )fe(k„ . . . )Ui{0)a^^{0, q+,\)\0) 



-fcj — fc- +ie 



mt. q states ^ 2q^ 

(0|a(k2, . . . )d(ki, . . . )Ui{0)d^{k,, . . . )|0) (0|d(k„ . . . )b{ko, . . . )Ui{0)a^^{0, q+,X)\0) 

int. q states 



E 



Q 



+(0|a(k2, . . . )d(ki, . . . )&(ko, . . . )Ui{0)a^\O, g+, A)|0) 



(8) 



where the second, third and fourth lines correspond respectively to the contributions of the diagrams ^a) , [2jb) and 
[2jc)-|-[2][d). For the diagram[2][a) (resp. Elb)), there is a summation over the momentum {kf , k^) and quantum number 
of the quark (resp. antiquark) on the internal line. Inserting the expressions (jA13[) . (jA14p . (jA17l) and (jAlSp for the 
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vertices and performing elementary algebra, one obtains the momentum space expression 



|7;(9^Q^A)^)^ 



e.9 



2(27r 



V^o Jo V^l JO ^2 



'lO,A2 

(1-zi) [l-2zi + 2feoA] (£A-ki) 
(1-^1-^2) (zi(l-zi)Q2+ki2) 
(1-zo) [l-2zo-2feoA] (£A-ko) 



d'kn 



d^ki /■ d^k: 



(27r)2 y (27r)2 J {2n) 



Z2 



1- 



(l-^l) 



l-2ho A2 



l + 2/io A2 



k2 ^ ki 

Z2 1-zi 

k2 , ko 



Z2 1-2^0 



JA,A2 



3A,-2ho 



2ho 



1-zo 



(27r)2,5(2)(ko+ki+k2) 



ZO Zl Z2 _ 

^ E ^/ E (^")aoAi^^(''"' ^0+' '^O' ^0' /) ^^(ki: kt, ~h,, A,, /) at(k2, fc2+, A2, a)|0) . (9) 

/ Ao,Ai,a 



In the case of a longitudinal photon, one needs to always start with the effective vertex (jA21l) . Since the latter 
already corresponds to a piece in a current-current instantaneous coulombian interaction, as explained in the appendix 
A 3[ there is no diagram with instantaneous splitting of a longitudinal photon into a quark-antiquark-gluon system. 
Hence, in the longitudinal photon case, one has only diagrams analog to^a) and[5Jb), so that 



6t(ko,...)d^(ki,...y(k2,...)|0) 



qqg states 



^0 ^1 ^2 



E 



|a(k2,...)6(ko,...)^/(0)&t(k»,...)|0) V^. 



-E 



mt. q states 

'"'a(k2, . . . )d(ki, . . . )Z^/(0)dt(ki, . . . )|0) V^2(,+ ,Q2)_,g(k„,... ) g(k.,... ) 



-") 



( 

int. q states 1 2g+ 

which gives, after substitution with the appropriate expressions, 
' ^ '"'wg ^ 7o Z2 J {2tt)^ J 27r 2 

'lO,A2 

* ( k2 I ki 

^A2 V 22 ^ 



(10) 



d^ko /■ d^ki f d^ka (2^)2^(2) (ko+ki+k2 



(2 



\^ ZO Zi Z2 



{l-ziVzi 

(1-Z1-Z2) 

(l-zo)2zo 



^^2 
(l-^l) 



l-2/io A2 



I + 2/10A2 



(zi(l-zi)Q2 + ki2) 









^A2- 




1-zo y 



(zo(l-2o)02 + ko2) 



(I-Z0-Z2) 

E^/ E (r'')^^^^6^(ko,fc+ /io,Ao,/)dt(ki,fc+,-/io,Ai,/)at(k2,fc2+,A2,a)|0). (11) 

The mixed space representation of those quark-antiquark-gluon components of the transverse and longitudinal 
photon wave- functions are obtained by Fourier transforming the creation operators as in (jAlip . In both cases, one 
gets an expression of the type 



^ dzo dzi dz2 , . 
—= / —= / 0(20 + 2^1 + 2:2-1) 



d xo /■ d xi /■ d X2 



(27r)2 J (27r)2 J {2n) 



eg 

2(27r) Jo Jq y/z{ Jq Z2 

E ^?!i(xo,Xi,X2,Zo,2l,Z2,/io,A2,(A)) 51^/ E (^'')aoA, 
ho,\2 f Ao,Ai,a 

X fe^(xo, k^, ho, Aq, /) d\xi,kf,-hQ, Ai,f) a''{x2,k^, A2, a)|0) , 



(12) 
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where, using the integrals (|B3I) . (|B4I) and (IBS 



$1?® (^Xo , Xi , X2 , , -21 , ^2 , ^0 , A2 , 



A = 



QXKi(gX) 



X2 



x<^ zi(l-zi) [l-2zi + 2/io A] £a- (^xio-Y7^X2o 



Z2 



-zo(l-zo) [l-2zo-2/io A] £a- xoi 



-Zo Zi Z2 (5a, A2 



Z2 



1-zo 



X21 



'J\,-2hQ "X,2ho 



1-Zi 1-Zo 



1- 



^^2 



(1-zi) V 2 

Z2 ( l + 2/io A2 

(1-20) 



£a,-x2o 



eA.-x2i 



1-21 



(13) 



and 



i>|''^(^Xo,Xi,X2,2:o,2l,Z2,^0, A2) = 2iQ Kq (QX) i Zi (1 - Zi) 



1- 



^^2 



(1-zl) 



-zo(l-zo) 



1- 



2:2 /l + 2/ioA2 



(1-20) 



l-2/ioA2 



eL-X2i 



''20 



In order to get more compact expressions, the variable X has been introduced, which is defined as 



X2 = Zi(l-Zi) Xio- 



Z2 



zo(l-zo) xo 



1-zi 

Z2 

1-zo 



X20 



X21 



22(1-21-^2) 

(1-21) 

22(1-20-^2) 
(l-^o) 



•-20 



= zi Zo a;io + ^2 2^20 + ^2 21 2:21 
Those three expressions are indeed equal when zq + zi + Z2 = 1 . 



(14) 



(15) 

(16) 
(17) 



Some comments about kinematical effects in the results 



Let us discuss the variable X and its various expressions in more details. The expression (1151) is the one naturally 
obtained for the argument of the Bessel function Ki(a;) (or Ko(a;) in the longitudinal case) when calculating the 
contribution of the diagrams [5](a) andt^fc) thanks to the integrals (jB3[) and (jB4p (or (jB5[) ). whereas the expression 
([TS]) is the one naturally obtained from the diagrams [2Kb) andj^Jd). It is quite remarkable that those two complicated 
expressions are actually equal when Z0 + Z1 + Z2 = 1. 



1. Recoil effects 

If one makes the approximation ki ~ — kp in the energy denominator ((5^ + kg/zo + k^/zi + k2/z2) of the momentum 
space expression (jO)), one gets the simplification 

•22 

xiQ-— — X20 xio (18) 

1 Z\ 

in the expression of X in the contribution of the diagram [^Ja) to the mixed space result. In momentum space, that 
approximation clearly corresponds to neglect some recoil effects and thus keep the same transverse momentum for 
the quark before and after the emission of the gluon. Hence, for the the diagram [^Ja), it is natural to introduce the 
transverse point xq' , defined by 

xio' = xio - 7^^x20 , (19) 

1-21 

and interpret it as the would be position of the quark in the absence of gluon emission. Since the point xi drops from 
that definition and 1 — zi = Z0 + Z2, the equation (|19p is equivalent to 



(20 + 22) Xo' = Zo Xo + Z2 X2 



(20) 
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Notice that the points (or Xq' ) are the positions (or would be position) of the partons at the time j:+ = when 
they interact with the target, but not at other values of x~^. 

On the other side, in light cone coordinates, the Poincare algebra has a transverse Galilean subalgebra where V" 
play the role of energy operator and thus of time, plays the role of mass and V± of momentum. This formalism 
then associates a transverse Galilean velocity v = k/A:"*" to a free massless particle of momentum and k. Then, 
the conservation of transverse momentum at the gluon emission vertex in the diagram [2l^a) implies the relation 

(■Z0 + Z2) Vq. = Zo Vq + Z2 V2 , (21) 

between the transverse Galilean velocity of the parent quark vq' and the ones of the daughter quark and gluon vq 
and V2, which is strongly reminiscent of the equation ([^0]). 

In our perturbative description of the projectile, the partons should be free streaming in the transverse plane 
between two successive interaction vertices, or between the last vertex and the interaction with the target at = 0. 
Hence, geometrically, the diagram [2][a) looks as in FigEl In the classical limit, there is a relation 

X, = x„ + (Ax+) V, , (22) 

for i 0, 2 or 0', where the emission of the gluon by the quark happens at a time — — Aa;+ and a transverse 
position x^. That relation (j22|) is sufficient to make the equations ([20| and ([21]) equivalent. Beyond the classical 
limit, a relation of the type (1221) should hold when using a basis of wave packets to diagonalize the free hamiltonian 
instead of the momentum space or of the mixed space Fock states. In that case, the formula (1^^ would relate the 
position of the center of the wave packet to the average momentum. 

Hence, both by tracking the origin of the expression (IT51) in our calculations and by discussing the Galilean subal- 
gebra of the Poincare algebra on the light front, one finds that the relation ([^0]) is the position space version of recoil 
effects associated to the gluon emission. This discussion of course generalizes to other diagrams and to the expression 

m- 

Generically, for any initial state parton cascade described in mixed space, we thus know how to take recoil effects 
into account at each splitting. Consider a parton which splits into two in this cascade, which has a longitudinal 
momentum fraction Zparont and which would be at the position Xparcnt at x"^ = if it was not splitting before. Then, 
the two daughter partons must have longitudinal momentum fractions Zdi and Zd2 and positions Xdi and Xd2 at 
x'^ — satisfying 

^parent = Zdl + Zd2 (23) 
^parent Xparcnt = 2dl Xdl + 2d2 Xd2 , (24) 

in order to conserve longitudinal and transverse momentum respectively. 
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2. Formation time of multi-parton fluctuations 



Thanks to that discussion of recoil effects, it is natural to rewrite the equations ((T5|) and ((T6| as 

X' = z,(l-z,)a:l,, + -^^xl, (25) 



{z2 + Zq) 
{Z2 + Zl) 



zo{l~zo)x\,^ + , s a^ai , (26) 



where xi' = {zi'x.i+ Z2'^2) I (zi + Z2) is the would be position of the antiquark at x'^ = in the absence of gluon 
emission in the diagram [2jb). In the expression ([25]), which is obtained for the diagram [SJa), the first term is 
obviously associated with the photon to dipole splitting and is a direct analog to the term appearing in the argument 
of the Bessel functions in the LO wave functions (O and (jG]), whereas the second term is associated with the emission 
of the gluon by the quark. And similarly, in the expression (1^51) obtained from the diagram [^l^b), the first term is 
again associated with the photon to dipole splitting and the second to the gluon emission by the antiquark. 

It has often proven useful for the physical interpretation of perturbative results to associate a formation time (see 
6.17. [l^l) to a parton splitting. The formation time is basically the time it takes before the two daughter particles 
are away enough from each other compared to their wavelength so that they lose their quantum coherence and can 
thus act independently from each other. For the splitting of a parton of momentum (fc^i + ^^2'^) ^'^^'^ daughter 
partons of momentum (fc^j^, k) and (fc^2' ^he formation time writes 

'n^orm — '7T+ 1 + \ \ 2 ' \ ' ) 



Up to a factor 2(j+ , each of the terms in the expressions ((25)) and ([26|) is then the mixed space version of the formation 
time associated to a parton splitting, with the inverse of the relative transverse momentum of the daughters replaced 
by the their transverse distance at x^ = 0. The quantity 2q~^ can thus be interpreted as the sum of the formation 
times associated with the two vertices of the diagram[2Ja) (resp. [DJb)) if one uses the expression ([^5)) (resp. ([25)) ). 

A further hint for this interpretation is the following. The wave functions ()13p and ()14p are exponentially suppressed 
at large QX due to the Bessel functions, so that effectively, one has the restriction 

^^^X'l^^-^^r,,, (28) 

where t^. is the lifetime of the virtual photon. This inequality means that the two successive splittings in the diagrams 
[5](a) and [D(b) should happen within the lifetime of the parent virtual photon in order to give a relevant contribution. 
In the LO wave functions ([5]) and ([6]), the same interpretation holds. The only difference is the presence of a single 
splitting vertex instead of two. 

At this point, we have understood the physics behind the expressions ([T5)) and ([TB)) . related to the diagrams [UJa) 
and [H^b) respectively, and also the argument of the Bessel functions in the LO wave functions. The first surprise is 
then that the expressions ()15p and ()16p are actually equal (provided recoil effects are not neglected), both reducing 
to the simple formula ([T7| . The second surprise is that even the instantaneous diagrams [2jc) and WiA) give Bessel 
functions with the same argument, despite the fact that the interpretation of 2q+ X^ as the sum of the formation 
times for two successive splittings does not make sense in those cases. One is thus led to the conclusion that 2g+ X"^ 
is the global formation time for a quark-antiquark-gluon system from a photon, and that it does not depend on the 
production mechanism but just on the final state of the cascade at x^ = 0. Notice that none of the three terms in 
the expression ()17p . taken individually, has an interpretation as a formation time for one splitting: the recoil effects 
are lacking in the first term, and the longitudinal momentum fraction of the parent is missing in the denominator of 
the second and third terms. The final expression ()17p for X^ is fully symmetric with respect to permutations of the 
particles. The formation time is indeed a purely kinematical quantity, independent of the nature of the particles and 
of dynamical properties of the cascade like the color fiow. 

When analysing the structure of the perturbative contributions to the quark, antiquark plus n gluons component of 
the photon wave function, one can see that our previous discussion generalizes. At any order in light-front perturbation 
theory, there should be a Ki(Q A'„) (resp. Ko(QX„)) Bessel function which factorizes from the n-partons component 
of the transverse (resp. longitudinal) photon wave function, with a common variable Xn for all diagrams leading to 
the same n-partons Fock state. And the mixed space expression for the formation time of a system of n partons 
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should write 

n-l 

Tform, n part = 2g+ X'i^ = 2g+ ^ Zj X^^ . (29) 

For the 2 and 3 partons cases, that formula (j29l) is indeed given by our perturbative calculations. In the appendix [Cl 
hints for the validity of the expression (|29|) for the 4 particles case are provided. 



D. Virtual corrections 



So far, we have calculated the quark-antiquark-gluon component but not the quark-antiquark component of the NLO 
virtual photon wave-functions ([7]). The direct calculation of the latter in light-front perturbation theory would involve 
loop diagrams and the renormalization of the theory, which is particularly cumbersome in this context. However, 
those NLO quark-antiquark components can be easily determined according to a prescription based on probability 
conservation (see e.g. 0, LO, the normalization of the wave-function ^ is non-trivial because we have 

dropped its colorless components, containing for example the photon itself or a lepton pair. Then, at NLO, QCD 
corrections are allowed, but not electroweak corrections. Hence, the global normalization of the colorfull sector of the 
photon wave-functions is the same at LO and at NLO with respect to QCD, so that 



LO 



(7T,L('z'+,Q^(A))^7fL(9+,Q^(A))^ 

= (7t,l(9'+,Q',(A))^7t,l(9+,Q',(A))^) + (7T,L(g'+, Q', (A))^ 7t.l(9+, Q', (A))^) 

00 \ / 00 000 \ / 



qq \ "I qq qqg 

From the relation (U), one gets the normalization of the wave- functions at LO 
;7t.l('z'+,Q',(A))^7t.l('Z+,Q',(A))^ 



LO 



LO 



(/lo), (A) 



(30) 



(31) 



And the normalization of the NLO quark-antiquark-gluon component is obtained from Ea. (|12p as 



7fL('z'+,Q',(A))^7;,L('Z+,Q',(A))^ 



X 1- 



ma 

2 



.9+ 



(27r)2 



d xq / d xi / dzi 



a I / — 2^ |$:^^^Xo,Xi,X2,l-Zi-Z2,Zi,Z2,/lo,A2,(A) 



(32) 



with the notation 



a= — as = y:-^ . 33 



The expression ((32)) has divergences for X2 — >■ xq, for X2 — >■ Xi and for Z2 — )■ 0, which should be regularized in some 
way. However, we will keep the regularization implicit here, for convenience. The first two divergences will disappear 
in the expressions for the virtual photon cross section, and we will come back in section UlI Bl on the Z2 — >■ divergence, 
associated with the low x factorization and evolution. 

The quark-antiquark component of the NLO virtual photon wave-functions can be parameterized in a similar way 
as the LO one (g]) 

|7f.(.^g^(A)),)^^_ ^ il^r^^^-^--^) '^>^^.(xo,x„zo,z„.o,(a)) 

xE^/E bH^o,k+,ho,AoJ)d^{^i,k+,-ho,Ao,f)\0). (34) 

/ ^0 
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This leads to an expression analog to (pT|) for the normalization of that component. 

Then, let us follow the standard prescription 0, [13] and assume that the probability conservation relation ([50]) 
actually valid at the level of the integrand, i.e. for each value of Xq, Xi and zi 



IS 



|$|^^(xo,Xi,l-zi,zi,(/io),(A)) |«'t!l(xo,xi,1-zi,zi,/io,(A) 







"~ y" X! |*t!l(xo,xi,X2,1-zi-Z2,2i,^2,/io,A2,(A) 



(35) 



III. PHOTON CROSS SECTIONS AND IMPACT FACTORS AT NLO 
A. From wave-functions to cross sections 

Thanks to the general formula (jA6[) . one writes the longitudinal photon cross section cr^ [-^] on a given classical 
Shockwave field at NLO accuracy as 



U+ ) 



NLO 



(2-)^(^-l) I 



9? ??9 



'H 



(36) 



Here the relation ([50)) has been used. The eikonal scattering operator Se does not change the partonic content of 
the projectile and has purely real matrix elements for colorless projectiles. That is why one can drop the real part 
operator Re and split the Se matrix element at NLO into the qq and qqg sector contributions. The transverse virtual 
photon cross section crj is defined to be unpolarized. Hence, the expression for cr^^ [-^] at NLO accuracy is analog to 
the expression ([55)1 for cr^ except that one has to average over the photon helicity A. 

Thanks to the relations (jA3l) and (|34|) the NLO matrix elements oi Se writes in the qq sector 



;7fL(<z'+,Q^(A))^ 5^,|7fi(<z+,Q^(A))^; 

= 2ns(^-l) ^f^Y.^) I d^xo / d^xi5oi[^] f Az, Y 1$??^ (xo, xi, 1-zi, zi, /.o, (A)) f , (37) 

\q / [ ^) J.- J J Jo 

with the notation 

S,M] = tr (C/[^](xo) f/n-4](xi)) (38) 
for the so-called dipole S- matrix. Because of the relation (pSj) . the equation (|57)) can be rewritten as 
;7T,L(9'+,Q',(A))^|^s|7fL('Z+,Q',(A))^) 



99 



99 



2t:S 



Ot^^''/ '^''""I ||<i>^o(xo,Xi,l-zi,zi,/.o,(A)) 



1-21 



~ f~^Yl kTi(x0,Xl,X2,l-Zi-Z2,2:i,2:2,/l0,A2,(A) 



A2 



(39) 
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Using the relations (jA3p and (|T2|) . the calculation of the NLO matrix elements of Se in the qqg sector gives 



7T,L('7'+,Q',(A))^^i=;7T.L(9+,Q',(A))^) _ =27r^ M^-l 



4 Oip,rri, 



g+ ; (27r)2 ^ / 



d xo / d xi 



t "^'^ /' " { E [-^] (^2)],, tr (C/[^](xo) T-^ t/n-4](xi) T") } 

|$^'^f,(^Xo,Xi,X2, l-Zl-Z2,2l,Z2, /lo, A2, (A) 



= 27r5( ^-1 



2-/Vc CHem 2 



(27r)2 



/ 



d Xo / d Xi / d^i a 



d'x2 /-^^"^ d^ 
27r in ^2 



5o2[^] 52i[^]-^5oi[^] 



^ |$|?£(^xo,xi,X2, l-2;i-2;2,2;i,z2,/io, A2, (A) 



where the second expression is obtained thanks to the identities 

[V[A]i^)],, = 2tr(t/[^](x)r"C/t[^](x)T^ 
1 „ „ 1 



and 



AB\- J CD ~ ^^A..dSb,C- ^SA^BSc,D■ 



J.a^ 



(40) 
(41) 

(42) 



Finally, inserting the results (jSlJ , (|39|) and (|40| into the equation (pS)) or its analog for the transverse photon cross 
section, one gets the photon cross sections for a given classical gluon Shockwave 



(2^) 



2 A^^f 

f 

l-zi 



d Xo / d xi / dz 



dSc2 
2n 



dz2 

Z2 



l-Soi[A] Ij.,^(xo,xi,l-zi,zi) 



Sai [A] - Sq2 [A] S21 [A] Xj,^^^(xo,xi,X2, l-2;i-z2, zi, Z2) 



with the impact factors defined as 

lf'°(xo,Xi, 2:0,^1) = i ^ |$|;°(xo,xi,zo,zi 

ho 

l|;°(xo,xi,zo,2i) = ^ ^ |^'t'^(xo,xi,zo, 2:1,^0 

ho, A 

^°(xo,Xi,X2,Zo,Zi,Z2) = ^ E |*l'^(xo,Xi,X2,zo,zi,Z2,/io, 



'J0.A2 



I-f^°(xo,Xi,X2, 20,^1,2:2) = ^ E |*T^(^0,Xl,X2, 



Zo,Zi,Z2,/lo, A, A; 



(43) 

(44) 
(45) 
(46) 
(47) 



Then, from the expressions ([5]) and ^ , one recovers of course the classic results of Nikolaev and Zakharov [l| for 
the LO impact factors for the dipole factorization 



Ij^ (xo,xi,zo,-2i) = 4(3 ZqZ^KqIQJ zqZixI^ 



l|°(xo,xi,zo,2i) = [zq + zf]zQZiQ^Kli Q 



ZqZiXqi 



(48) 



(49) 



The NLO impact factors, which constitute the main results of the present work, are calculated using the intermediate 
results ([T^ and p^ . and write 



Z2 
1-Zi 



Zj^^O(xo,Xi,X2, 20,^1,^2) = 4Q^Kl{QX){zl{l~z,)^ +^o(l-^o) 



Z2 
l-Zf) 



''20 



1-21 



-2zi(l-zi)zo(l-zo) 



Z2 



Z2 



2{l-zi) 2(l-zo) 



X20 ■ X21 



(50) 
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and 



(Xo , Xi , X2 , 2;o , 21 , 22 ) 



QXKiiQX) 



zf(l-zi)^[z? + (l-zi)2 



Z2 



1-Zl 



X20 



+Zo2(l-zof [2o' + (l-2o)2] fxoi--^X2iy ^ 



l-zo J 



+2zi{l- zi)zo{l- zo)[zi(l- zq) + zo(l-zi)] 

ZO 21 ^1 (20-21)^ (x2oAX2l)^ 2 

' + Zo Zi Z2 



22 



22 



(l-zi)(l-zo) 



^2 ^2 



2(l-zi) 2(l-zo) 

ZqZi (1-2i)^ 



xio- 



22 



1-Zl 



(1-zi) (1-zo) 



22 \ / X20 

Xio--; ^20 M — 



+ Zq Zl Z2 



Zo Zl 



(l-2o)2 



(1-zo) (1-zi) 



22 \ / X21 

Xoi--; X21 • 

1-zo J \x^^ 



v2 ^2 ^2 



1-Zl 
1 



L(i-zi)2 (i^zoy 



2 75 



X20 ■ Xqi 



'-20 



22 

1-zo 



X21 



X20 • X21 



Xnri Xn 



(51) 



where x A y = e*-' x' y-' and 



Viz) 



l + (l-z)^ 



(52) 



The two first terms m the NLO mipact factors (|50p and ([STt come from the squares of the diagram Ufa) and (Hb) 
respectively (or their analog for longitudinal photon). Since V{z)/z is the quark to gluon DGLAP sphtting function, 
up to some color factor, those two terms correctly reproduce the first step in the LO DGLAP evolution of the virtual 
photon in the appropriate limits x|g — and — 0, relevant in the so-called resolved photon kinematics. By 
contrast, the first contribution to the more familiar DGLAP evolution of the target in DIS is not clearly apparent in 
the results (ISUl) and (ISTI) . although it should be somehow present there. In the limits x^q — and 0, the 

two first terms in the NLO impact factors (|50p and ((5T|) have the singular behavior responsible for the logarithmic 
divergences already mentioned in the section III Dl However, those divergences are regulated in the cross sections 
43]) by the property of the color transparency 



. [A 



for 



(53) 



Hence, the integration over X2 in the formula (j43[) does not need any regularization anymore. 



The third term in the NLO impact factors ([50]) and (|5T|) comes from the interference between the diagrams [2ja) 
andt^fb). However, in the transverse case (j5ip . the fourth term is an another contribution arising from the same 
interference. In intermediate stages of its calculation, this additional interference contribution is proportional to 
(xio' A Xi/o)(x2o A X21), where, following the notations of the section ITTCl xio' (resp. Xi/q) is the parent dipole vector 
in the diagram [2](a) (resp. [2](b)) including recoil effects. Notice that this contribution is the product of two terms 
which are odd by exchange of Xi and zi with Xq and zq, so that it is globally even, as are other contributions to the 
impact factors (l48l) . p9|) . (|50p and (|5T|) . Such an odd-odd contribution requires two vertices with a non-trivial tensor 
structure. In particular, the effective vertex (jA2ip or the instantaneous vertex (jAlSP are essentially scalar so that 
they cannot provide an odd product. All this explains why the odd-odd contributions start to appear at NLO for the 
transverse photon case, and not before NNLO for the longitudinal case. In addition, the antisymmetry of each factor 
implies that only interference terms can produce a non-vanishing odd-odd contribution. Recoil effects are also crucial 
for the appearance of the odd-odd contribution in the result ([5T|) . by ensuring that Xio' A xi'o Xio A Xio = 0. 

Finally, the last term in the expression ((5T|) comes from the square of the instantaneous diagrams [2l[c) and[2ljd) and 
the fifth and sixth terms from interferences between instantaneous and non-instantaneous diagrams, which explains 
why they have no analog in the longitudinal photon case (f50| . 



B. Subtracting high energy leading logs 

In order to get the final result for the photon-target cross sections from (|43|) . one still has to take the CGC average 
over the classical gluon field of the target and to deal with the logarithmic divergence at Z2 — >■ 0. For generic field 
A~ of the target, the dipole S-matrix suffers from the so called rapidity divergence. Those two divergence are 
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actually related and can be simultaneously regulated by a single cut-ofF preventing double counting of gluons. The 
idea is roughly to consider, in some particular frame, the right moving gluons as part of the projectile and the left 
moving gluons as part of the target. In practice, various choices of factorization scheme are possible, but the most 
convenient in our case is to use a cut-ofF kj' in fc+, following e.g. 3J|. The CGC average is then restricted to fields 

A~ which have no modes of larger than fc , and in the formula for the photon-target cross section the lower bound 
for the integration over Z2 is taken to he z f — j . 

The lowest possible value kl^^^ for the cut-off fc^ is set by the target. Indeed, from the point of view of the target, 
a radiated gluon should be a fluctuation of lifetime Xjk'^ shorter than the lifetime of the typical non-perturbative 
constituents of the target. Assuming that each of those constituents carry a sizable fraction of the momentum P~ of 
the target^ and are associated to a typical transverse scale Qo, that bound writes 

1 1 2P- , fc+ Ql Ql 

-that - = ^>-n.i„ = ^^=x^. (54) 

For a hadron target Qo should be a non-perturbative QCD scale, and for a large nucleus target Qo can be identified 
with the initial saturation scale (5^(0). 

From the formula p3|) . it is obvious that the cut-ofF Zf should be smaller than 1 — zi. Hence, Zf satisfies 

^^<^/^l-^i- (55) 

It is convenient to introduce the variable 

analogous to but distinct from a rapidity range, in order to specify the range left for the evolution of the target. That 
variable is bounded by 

The CGC average over the classical gluon field, associated to the cut-off A:t, is then noted (. . . 

All in all, one can deduce from the equation (|43|) the final result for the transverse and longitudinal virtual photon 
cross sections 



^ ' f 



y"d^xo y^l^^i J dzi |l|;'2(xo,xi,l-zi,zi) 



1 - (5oi)y+ + alog ( — —] [ (5oi-5o252i)y+ 

(-1-21 



( -ir^ {Sqi- So2S2i)y+ I — A2:^^°(xo,xi,x2, 21,2:2) > , (58) 
J 27r J^^ Z2 ' J 

where the high energy leading logs have been separated from the rest of the NLO corrections by defining 

AI^f;'^(xo,Xi,X2, 21,22) = ^^^^(Xo, Xi,X2, 1-21-22, 21, 22) -I^f;°(Xo,Xi,X2, 1-21, 21,0) (59) 

and noticing that 

I^£0(xo,Xl,X2, 1-21, 21,0) =- Xj:g(xo,Xl,l-2l,2l) ^ ^\ . (60) 

^20 ^21 



^ In the case of a target nucleus, P is defined as the momentum per nucleon rather than the total one. 
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The general structure of the NLO dipole factorization formula (|58p look similar to the one of the results of Balitsky 
and Chirilli for the high energy operator product expansion of two currents [361] . However, a more detailed comparison 
would require to perform a rather cumbersome Fourier transform of their result with respect to the positions of the 
two currents and the separation of the transverse and longitudinal photon contributions. 

As expected, the evolution with respect to a change of the factorization scale Zf (or equivalently is governed 
by the B-JIMWLK equations. Indeed, the B-JIMWLK evolution of (5oi)y+ writes 

9y+ (5oi)y+ = a. /^TZ^-r^V (5o252i-5oi)v+ , (61) 

J 02 2X 

and (502521 — iSoi)y+ satisfies a similar but more complicated equation. Hence, the derivative with respect to Yf of the 

bracket in the second line of the equation ((55)) vanishes up to terms of order 0{o?\ The evolution of (5oi — Sq2 S2\)y^ 

in the third line of the equation (|58p also brings only irrelevant NNLO contributions. Finally, one might replace the 
lower bound of the Z2-iiitegral of the third line by 0, since the integrand is now regular. 

In order to optimize the resummation of leading logs, one should obviously choose a factorization scale Zj close 
enough to the upper part of its allowed range (|55|) . Accordingly, one should evolve (501)^+ and (502521 —5oi)y+ with 

the B-JIMWLK equations over an interval r+ ~ log ((1-zi) QVa^Qo)- 

As a remark, notice that the dipole S-matrix (5oi)y+ used in the factorization formula (|58p or its LO approximation 
then depends on the Bjorken x or only through the combination W"^ = Q'^ /x, as advocated e.g. in Refs.[45l. l46j. 

IV. CONCLUSION 

The main result presented here is the NLO generalization (|58p of the dipole factorization formula for DIS structure 
functions at low x, with the longitudinal and transverse NLO impact factors (|5n|) and (|5T|) . In ref. [s^, results 
in principle equivalent have been already presented. However, in that paper, the results were given only in full 
position space for the incoming and outgoing photons, so that they were not ready for use in phenomenology before 
cumbersome Fourier transforms, by contrast to the results presented here. In future studies, it would be useful to 
perform those Fourier transforms, so that the independent NLO calculations of ref. [s^ and of the present paper 
would cross-check each other. 

Obviously, the results presented here pave the way towards phenomenological studies in the spirit of refs. (3ll - [33| . 
but at NLO accuracy. For consistency, however, one might have to use the present results together with e.g. the NLL 
BK equation [s^ . |35| , whose numerical simulation seems to be a formidable task. 

The calculations have been performed under the assumption that the target is dense, and thus include gluon 
saturation effects. However, in the 2 gluons exchange approximation, they should provide the appropriate NLO 
factorization formula to be used with the BFKL equation in mixed space. 

As an intermediate step, the quark-antiquark-gluon components of the light front wave functions of virtual photons 
have been obtained, see equations (IT2]) . (IT3l) and (IT4l) . Those should be one of the building blocks for the NLO dipole 
factorization of other DIS observables, typically less inclusive. 

With DIS phenomenology in mind, it would be useful to extend the results of the present study to quark mass 
effects. Indeed, the charm quark is known to provide a sizable contribution to DIS structure functions, and the bottom 
quark a non- negligible one. Quark mass effects are also required for phenomenology at low . 

As a byproduct, the results presented here provide an interesting insight into the kinematics of initial state parton 
cascades in mixed space. It has been inferred that the prescription ([M]) allows to take recoil effects exactly into 
account in mixed space. Finally, the generic formula (1291) for the mixed space expression for the formation time of a 
arbitrary n-partons state has been conjectured. That expression depends only on the final state of the shower, not 
on the precise diagram leading to that final state. 
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Appendix A: Light-front perturbation theory for DIS at low x 
1. Eikonal scattering and total cross section of a dilute projectile on a Color Glass Condensate 



Let us consider a right-moving relativistic projectile scattering on a very dense and highly boosted left-moving 
target. In DIS at low x and similar processes, only the gluons with low momentum fraction of the target are probed. 
For a dense enough target or at low enough x, the occupation number of those soft gluons becomes non-perturbatively 
large, so that the state of the target is more appropriately described in terms of semi-classical gluon fields than of Fock 
states with a finite number of gluons. This is the starting point [23l - l25l | of the Color Glass Condensate (CGC) effective 
theory for the target and of related frameworks jl^]. For the collision of a dilute projectile on such a dense target, the 
disturbance of the semi-classical field of the target by the interaction is a higher order effect, so that only diagonal 
terms in the density matrix for the soft gluon field of the target are probed and thus the quantum average over the 
target gluon field reduces to a classical statistical average in the CGC formalism, up to perturbatively calculable 
corrections. 

Then, one has to calculate the high-energy scattering of the projectile on an arbitrary classical gluon field. This 
can be done with the light-front formalism of ref.[39| generalized to QCD. In the limit of infinite boost, the target 
field-strength is Lorentz contracted to a Shockwave localized at'^ x'^ = 0. Hence, the interaction of the projectile with 
the classical field happens during an infinitely short x~^ time interval around x'^ ~ 0. In the Heisenberg picture, the 
projectile is thus in a state \iH) for x"*" < and switches to a state S'^jiij) for x'^ > 0, where Se is a unitary operator 
which depends on the classical gluon field. 

In the high-energy limit, the typical time-scale along a;"*" for interactions between the partons in the projectile 
becomes infinitely larger than the x"*" time interval taken by the projectile to pass through the classical field, suggesting 
that these two types of interactions factorize from each other. Thus, different partons present in the projectile at 
a;"*" = should scatter independently from each other on the Shockwave field. The transverse position x of partons 
stays obviously constant during their infinitely fast interaction with the infinitely thin Shockwave, but their transverse 
momentum k is typically modified. Moreover, in the limit of infinite boost of the target, the Shockwave field is 
independent of x~ , so that it contains only gluon modes with fc+ — 0. Hence, the momentum of partons from 
the projectile cannot be modified by scattering on the Shockwave field. It is thus convenient to use a mixed space 
representation, specifying the longitudinal momentum and the transverse position x of partons in the projectile. 
We then use a Fock state basis to describe the state of the projectile at a;"*" = 0, just before or just after the collision. 
The normalization of the annihilation operators a, b and d for gluons, quarks and anti-quarks is chosen such that 



{6(x',fc'+,/i',A',/'),6t(x, 



,/)} = [d{^',k'+,h',A',ndH^,k+,h,Aj)} 

= (2^)32fc+ S{k+-k'+) (5(2)(x-x') S'^^'^' 5^^^' 5f'f' 
a(x', fc'+. A', a'), at(x, fc+. A, a)] = {2TTf2k+ 5{k+-k'^) (5(2) (x-x') 5^'^' S""'"' , 



(Al) 
(A2) 



where h — ±1/2 and A = ±1 are the helicities, A and a the fundamental and adjoint color indices, and / the fiavor. 
At high energy, the scattering of partons off the gluon Shockwave is eikonal, so that Se acts on Fock states by only 
color rotating each partons by a Wilson line defined along its trajectory through the Shockwave, i.e. 

Se b\-Ki,k^,hi,Ai, fi)--- d^(x™, /i™, A™, /„) • • • a1'(x„, fc+, /i„, A„, /„)••• |0) 
= E E E [-4] ■■■[U^ [A (X„0] A^B^ ■■■[V [A] (X„)],„,„ • • • 



Al,... A^,... a„,... 

xb\xi,k+,hi,AiJi)---d^:ic. 



^ ' : -^^L : J L J ' ' ' y.^m ; ; ^mi 5 fm) ' ' ' O,^ (^n : k^ , A^ , fln) ' ' ' 1 0) ; 

where the fundamental and adjoint Wilson lines are respectively defined as the path-ordered exponential 



U [A] (x) = V exp 
V [A] (x) = V exp 



ig J dx+ T°y^-(a;+,x,0) 
ig J dx+ Aaix+,x,0) 



(A3) 

(A4) 
(A5) 



In this study, the hght-cone coordinates are defined as = (x" it x^)/\/2. 
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involving the classical Shockwave gluon field A'^ {x~^ , x, a; ) . 

Finally, the projectile-shockwave total cross-section is obtained by the optical theorem, which writes in this case 

Re({tH{q'^)\l-SE\tH{q+))) 

\q^i2.)Siq'--q^) ' ^^'^ 

where q+ is the momentum of the projectile. 



2. Light-front wave-function of the projectile in perturbation theory 

In order to calculate the total cross section of the projectile on a given Shockwave field using ea. (IA6|) . we now only 
need to know how to write the incoming Heisenberg state of the projectile in terms of the Fock state basis for free 
on-shell partons at x+ ~ 0. This can be done perturbatively in the framework of light- front wave functions [39| . The 
light-front hamiltonian writes =T + U, where T is the free part and li the interaction part (see ref.ji^ for 
more details). In the interaction picture, the evolution of operators is generated by T, for example 

lli{x+) = e''^^^ Ui{0) e-*"^^^ , (A7) 
for the interaction operator, and the states \ii{x'^)) evolve as 

\ii{x^)) expi~i [ \x+Ui{x+)] \ii{x+)). (A8) 



The interaction picture is defined here in such a way that it coincides at x'^ = with the Heisenberg picture, i.e. 
|z/(0)) = Then, expanding the exponential in (jA8[) for X2 — and x^ — > ~oo and inserting several times the 

decomposition of the identity 



J" 

over a basis of Fock states \J-) at a;+ ~ 0, one finds 



1 



n=l J^„ J^i 



where Tjr is the eigenvalue of the free hamiltonian T corresponding to the state \ J^). Hence, we are forced to choose 
the standard momentum-space Fock basis instead of the mixed-space Fock basis considered previously, which does 
not diagonalize 7". The two Fock bases are simply related by transverse Fourier transform of all the creation operators 
present in the state, with 

6t(k, h,AJ)=J^ e'''- b\^, h, A, f) (All) 

and similar relations for S and . 

Thus, the general method for the calculation is the following. First, calculate the momentum-space light-front wave 
function of the projectile {J^Iih) to the appropriate order in perturbation theory using the general expression (jAlOp . 
Then, perform the required Fourier transforms in order to get the light-front wave function in mixed-space. And 
finally, use the relations (jA3[) and (|A6I) or analog ones, in order to obtain the total cross section or other observables 
for the scattering of the projectile off a given classical gluon field, before performing the CGC statistical average over 
the gluon field of the target. 

Let us come back to the formula (jAlOp in order to specify the last details. First, the eigenvalue Tjr is simply the 
sum over the k~ for the partons i present in the Fock state. Remember that in light-front perturbation theory, only 
physical on-shell partons are included in the Hilbert space, so that k~ = ki^/(2fc^) for massless partons. 
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In the formal sum over Fock states used in equations (jA9l) and (jA10[) . there is a summation over the number of 
partons of each type present in the Fock state. And for each parton present, there is a summation over its quantum 
numbers and a phase-space integration 



-°° dk+ 



(27r)(2fc+) J (2^) 



d^k 



(A12) 



consistent with the normahzation of creation operators. 

Inserting the expressions for the quantized free fields at a;"*" = into the interaction part of the hamiltonian Ui{0), 
one can obtain the required expressions for the vertices. In this study, we need two types of QCD vertices on the 
light-front: the quark to quark and gluon splitting and the antiquark to antiquark and gluon splitting, which write 
respectively 



i{k",k"^,X,a)b{k',k'^,h',A',f)Ui{0)b\k,k+,h,AJ)\0) = (27r)35(fc'+ + fc"+-fc+) 5(2) (k'+k"-k) 



el ■ 



k" / l + (2fe)A \ k^ / l-(2/i)A \ k^ 



J k+ 



(A13) 



and 



i(k", A, a)d(k', h\ A', f')Ui{0)d^k, k+,h, A, /)|0) = (27r)3j(fc'+ + fc"+-fc+) J^^) (k' + k"-k) 



k'- 



A k^ n-{2h)X \ k_ 
J k^'\ 2 J k^ 



(AM) 



in the case of massless quarks. Here, ex is the transverse polarization vector for transverse gauge bosons of helicity 
A, i.e. 



£x = 



1 / 1 



V2 



(A15) 



which satisfies the relations 



^ eX ei = 5'^ and 



A=±l 



A=±l 



^A — * 



(A16) 



and e^^ is antisymmetric with e^^ = -1-1. 

We also need the QED vertex for the splitting of a photon into a quark and antiquark dipole 

(0|d(ki, hi,Ai,fi)b{ko, k+,ho, Ao, fo)Ui{0)a^HK k+,X)\0) = (2^)3^(fco+ + fci+-fc+) 6^^^ (kg-fki-k) 

A A A \^ f l + {2h„)X \ ki f l-{2ho)X \ kp 
X ee/o dfgj, dAo,Ai ^ho.-h, y4:ko k^ ex ■ — - j ^ j -j-f 



(A17) 



and the mixed QED/QCD vertex for the instantaneous splitting of a photon into a quark, an antiquark and a gluon 
(0|a(k2, fc+, A2, a)d(ki, fc+, /ii, Ai, /i)6(ko, fc+, ho, Aq, fo)Ui{0)a^\k, A)|0) = (2^)35(fco+ + fci + + fc2+-fc+) 



(5^2^ (ko + ki-|-k2-k) ee/o Sj^j^ (^°)aoAi '^''0 -/h ^a,Ao A/4fc^fc 



'^1 



(A18) 



where e is the proton electric charge and e/ the fractional charge of the quark of flavor /. 



3. DIS case: from lepton to virtual photon scattering 



In DIS, the initial asymptotic state of the projectile 00)) contains only one lepton, and |J^o) = {iii^oo)). As 
usual, only the LO contribution to DIS with respect to QED is relevant: the initial lepton couples to quarks via 
a single photon exchange. To that order, one can trivially factor out the leptonic part of the diagram in covariant 
perturbation theory, and then one has to study the scattering off the target of the intermediate virtual photon which 
can be transverse or longitudinal. In our formalism based on light-front perturbation theory, there are neither virtual 
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particles nor longitudinal photons in the Hilbert space, so that the factorization of the leptonic tensor is much less 
obvious but still holds, for the following reasons. 

In all energy denominators of the formula (jA10|) . the momentum (/cj^,k/) of the incoming lepton appears via 
Tjfq = = kf/{2k^). And the scattered lepton is present in all of the Fock states corresponding to intermediate 
or final steps in (|A10|) . in addition to the emitted photon or to colored partons. Hence the second term Tjr^ (with 
TO > 0) in all the energy denominators contains the contribution fcj" = kj /{2k'i'^) of the scattered lepton. Because of 
momentum conservation, the momentum of the emitted photon is = k'^ — k'i'^ and q = k/— kj. Although it cannot 
be interpreted as photon virtuality in this light-front formalism, it is still convenient to introduce the variable 

Q' = -{k^^kr){k^-k[^). (A19) 

In a frame in which q = 0, the contribution of the incoming and scattered leptons to any energy denominator writes 

k--kr^-^. (A20) 

Apart from that, the incoming and scattered leptons are involved in the formula (jA10| only in the first interaction 
vertex. Actually, DIS processes can be initiated in two ways in the present formalism: either the lepton emits a 
transverse photon, which then propagates and splits later into a quark antiquark dipole, or the quark antiquark pair 
production occurs in one step via instantaneous coulombian interaction between the quark and leptonic currents. The 
latter case can be understood as involving an instantaneous longitudinal photon exchange, so that this photon is not 
present in any of the intermediate Fock states J-m- Thanks to the integration over the azimuthal angle of the outgoing 
lepton performed for DIS, the interference between those transverse and longitudinal contributions vanish, so that 
they add up incoherently in the DIS cross section. 

The leptonic part of the contribution with transverse photon can be factorized out easily: just remove the first 
step in (|A10I) . including the lepton to lepton and photon vertex and the energy denominator for the lepton-photon 
intermediate Fock state. Then, one can use all the formalism presented before, but using the transverse photon instead 
of the lepton as projectile and replacing Tf^^ by — Q^/(2(7+) in (jAlOp . and calculate the total cross section crj for the 
scattering of the transverse photon off the target. 

By contrast, for the longitudinal contribution, one has to separate by hand the lepton to lepton quark and anti- 
quark vertex into a factor associated to the QED coupling of the leptons, a factor analog to an energy denominator, 
and a factor associated to the QED pair production of quarks. This last factor is interpreted as a fictitious vertex for 
the splitting of a longitudinal virtual photon into a quark anti-quark dipole, and writes 

V^.(,+ ^Q2)^,(k,fe+^^,^,/),-(k',fe'+,v,A',/') = (2iTf5{k++k'+^q+)6^'H^+'^') Sfj'dh'^-hSA,A'eef y/4k+k'+ ^ . (A21) 

Then, the total cross section a'l for the scattering of the longitudinal photon off the target is calculated using this 
splitting vertex instead of the transverse one (|JT7)) and the replacement of Tf^ by -Q'^/(2q+) in energy denominators. 

Hence, the hadronic part of the diagrams relevant to DIS can be calculated separately from the leptonic part in 
light-front perturbation theory as well as in covariant perturbation theory, leading to the photon target cross sections 
a'^ and cr^. And the full DIS cross section in the one photon exchange approximation is obtained via the same 
expression as in covariant perturbation theory 

w ^ ^ { h4) "?(--«^' + "J(^-<?^'} ^ <*^^) 

The variable y is defined as 

y = —. (A23) 

xs 

s is the Mandelstam s variable for the lepton-target collision, and the Bjorken x variable can be approximated for our 
purposes as 

2{q-P)- 2q+P-- ^^^^> 

The energy-momentum momentum of the target has a big component P~ , whereas P = by choice of frame, and 
P+ is given by the on-shellness condition = P^P^ = 2P~ P+. 
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Naively, cr'^{x,Q^) and aKxjQ"^) look independant of the Bjorken x in our formalism. Indeed, jij^) describes the 
projectile independently of the target, the eikonal scattering operator Se seems independent of the collision energy, 
as long as it is large enough for the formalism to be valid, and the CGC average over the gluon field of the target 
seems independent of the projectile. However, the contributions to {in) containing gluons at = have the usual 
divergence of Bremsstrahlung in the soft gluon limit. And the Wilson lines have rapidity divergences for generic field of 
the target. These two types divergences can be regulated by a common cut-off in longitudinal momentum or rapidity, 
as explained in section IIII Bl Then, one has a renormalization group evolution when moving this cut-off, which is the 
B-JIMWLK evolution in general, and reduces to the BK equation in a mean field approximation and to the BFKL 
equation in the dilute target limit. Then a^{x,Q'^) and a1{x,Q'^) become dependent on x through the factorization 
scale associated to that evolution. 



Appendix B: A few integrals 

When calculating the wave- functions for transverse and longitudinal virtual photons in mixed space up to NLO, 
one encounters the following integrals: 



j2. ikx j _ _ 

— =5 -k^ - z-QKi(Qlxl) (Bl) 

27r n j_v2 X 



-f k2 
d^k e*''-^ 



2^ Q%k2 

d^ki fd'k, e^--->+'^^--^« kj (^ + ^) 

(z,(i-zOQ^+k,^) (o^+fas+§+g 



Ko(Q|x|) (B2) 



2tt / 2tt 



(^^i.-j-^^^io) ^ J, (B3) 

d'ki fd^ka e»ki.xio+*k2X2o OXKiiQX) , ^, 

- Z2Zi[l-Zi-Z2) Tp^ (B4) 



27r ./ 2Tr f ^9. , (ki+k.)^ , kj ^ kjA "^-^^^ 

' (1— Zi— Z2) 21 Z2 



]2, 



1^1 i2i piki-xio-|-ik2 X2o (Sj_JSi1i - „ 



Ko(QX) . (B5) 



'^20 



Here, X is defined by Eq. pS)) . All those integrals can be performed by using the Schwinger representation of the 
denominators in order to transform the integrals over the transverse momentums into gaussian ones. 



Appendix C: Formation time for 4-particIe states 



Excluding diagrams with instantaneous interactions, there are two generic topologies for tree level diagrams de- 
scribing the splitting of one particle into four, up to permutations of particles. Those two topologies are illustrated 
in FiglH Following the prescription (p4|) for recoil effects, one writes in the case of the topology (A) 

Z2X2-HZ3X3 .„^. 

X2' = (CI) 

2:2 + Z3 

ZoXq + (Z2 + Z3)X2' 20X0 + ^2X2+2:3X3 

Xo' - ; ; - ; ; . {(ui) 

Z0 + Z2 + Z2 Za + Z2 + Z2 
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X-+ = 




FIG. 4: The two possible topologies (A) and (B) for diagrams without instantaneous interactions leading to a 4-partons Fock 
state. 



Then, one constructs the variable j-^j by summing the formation times in mixed space associated with each of the 
three vertices 



X 



4,M) = zi (zo + Z2 + Z3)a;o'i + \ \ 2^02' + — \ — 2:23 

^0 + ^2 + ^3 -,- 



22 + 2:3 

, ZO Xo + 22 X2 + 2:3 X3 ■ ^ 

21 (20+22 + 23) I ^ ■ Xi 



20 + 22 + 23 

In an analogous way, one has for the topology (B 

Xo' 
Xl' 



20 + 22 + 22 
20 (22 + 23 ) / 22 X2 + 23 X3 

Xo - 



22+23 



20 Xo + 22 X2 

20 + 22 

Zi Xl + Z3 X3 

21 + 23 



22 23 2 
22 + 23 



And the variable A"^ sums the formation time for the three splitting present for that topology 

v2 I t \( t \ "2 , Z0Z2 2 , Z1Z3 2 
^A,{B) = (2o + 22)(2i+Z3) Xo'i/ H XnQ H \ X 



= (2o + 22)(2i + Z3) 



09 T" ^13 
20 + 22 21+^3 " 



Zo Xo + 22 X2 21 Xl + Z3 X3 



20 + 22 



21+23 



20 22 2 



20+22 21+Z3 " 



2l 23 2 



(C3) 

(C4) 
(C5) 



(C6) 



It is then possible, using the constraint Z0 + 21 + Z2 + 23 = 1, to show that both X^ j.^^ and X^ f^g^ reduce to the 



expression X| given by the equation ([^ 

Xj (^N = Xl tg-. = Xj ^ Zq Zl Xoi + 2o 22 XI2 + 2o Z3 Xo3 + 2l Z2 x\2 + 2l Z3 + Z2 23 X^^ . 



(C7) 



Hence, even for larger n, the variable X^ associated to the formation time of a n-partons state seems still independent 
on the path taken by the parton cascade. The diagrams with instantaneous interactions should not destroy that 
property, like in the 3 partons case. Indeed, due to gauge invariance, they should not feature so different kinematic 
properties than non-instantaneous diagrams. 
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